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Abstract: We explore a mechanism for radiatively generating neutrino Majorana masses 
in a 5 dimensional orbifold SU{5) unification model without introducing right-handed 
singlets. The model is non-super symmetric and the extra dimension is compactified via a 
Si/{Z2xZ2) orbifold geometry. The necessary lepton number violating interaction arises 
from the Yukawa interactions either between a 10-plet or a 15-plet bulk scalar field and 
the fermion quintuplets which are residents on the SU (5) symmetrical brane located at 
one of the orbifold fixed points. The model is engineered to give realistic charged fermion 
masses and mixing and in the same time avoiding the rapid proton and neutron decays by 
geometric construction. The gauge unification can be maintained by adding extra fermion 
or scalar fields. The unification scale is found to be larger then 10^^ GeV by adding a 
bulk vector decuplet pair whose zero mode has masses around 10 ~ 100 TeV range. We 
found that neutrino mass matrix of the normal hierarchy type is favored by using 15-plet 
scalar. We give a solution of this type which has detectable /i — > 3e transition. On the 
other hand, by introducing 10-plet scalar, the leading neutrino mass matrix can only be 
inverted hierarchical and gives at most bi-maximal mixing. 
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1. Introduction 



Recent measurements of atmospheric and solar neutrino fluxes at the Super Kamiokande Q 
and the Sudbury Neutrino Observatory ||2| have provided compelling evidence for neutrino 
masses and neutrino oscillations. This received further support from the reactor KamLand 
experiment Q. Furthermore, the Wilkinson Microwave Anisotropy Probe Q imposes the 
constraint that the sum of neutrino masses to be less than .75 eV. Such a small value 
for neutrino masses is generally considered to be a harbinger of new physics beyond the 
standard model (SM) and the existence of a new scale between the Fermi and the Planck 
scale. In particular if the three neutrinos involved in weak interactions are Majorana in 
nature then clearly new physics is at play. The most popular suggestion of generating 
neutrino masses in the milli-electronvolt range is grand unified theories (GUTs) via the 
seesaw mechanism with or without super symmetry. Central to this idea is the introduction 
of one right-handed singlet neutrino per family of the SM fermions with a mass near the 
GUT scale. This is natural in SO{W) models since its fundamental 16 representation 
encompasses this singlet with the 15 fermions of the SM. For a recent review of neutrino 
masses in grand unified models see [Q. On the other hand small Dirac neutrino masses 
is considered unnatural due to the extreme fine tuning required. However, in theories 
with extra dimensions this can be generated by allowing the singlet neutrinos to be bulk 
fields. A small Yukawa coupling can be obtained due to the volume dilution factor if the 
extra dimensions are sufficiently large |^]. In both cases right-handed singlet fields Nji are 
necessary. 

In this paper we study the construction of neutrino mass without the benefit of 
in the context of grand unified SU{5) models with the minimal particle content^. This is 
a fundamentally different mechanism from the above mentioned constructions. Since the 
neutrinos in this scenario are Weyl particles of the SM and the resulting neutrino mass 
matrix is necessarily Majorana. In conventional four dimensional (4D) field theories one 
can use SM Higgs triplet to achieve this. However, the vacuum expectation value (VEV) 
of the triplet must be fine tuned to small values in order not to upset the highly successful 
custodial SU{2) relation of the SM gauge bosons as well as to generate of a sufficiently 
small neutrino mass. A second method is to radiatively generate neutrino masses at 1-loop. 
The prototype model was constructed some time ago |^ and the crucial ingredient is the 
introduction of a SU (2) singlet scalar field with non trivial weak hypercharge. The original 
version of the model gives a 3 x 3 neutrino mass matrix with zero diagonal elements and thus 
leading to bi-maximal neutrino mixings [^. This is ruled out by the data. However, simple 
phenomenological modifications can bring it to agree with observations |l^. All these 
constructions suffer from being rather ad hoc. It will be interesting if one can incorporate 
these attempts into the theoretically well motivated unification models. A more modern 
formulation of this makes use of progress in recent works on extra dimensions and the 
brane world scenario. It has been shown that the technique of orbifold projections applied 

""^It is also possible to generate neutrino masses without using Nr by R-parity violating interaction in 
minimal supersymmetric standard model|Q. Here we concentrate on non-supersymmetric models and leave 
questions such as quantum stability of scalar fields and naturalness issues as unsolved. 
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to GUTs [ 11 1 can solve some of the long standing problems such as doublet-triplet splitting 
that plagued 4D SU{5). This is further applied to the flavor problem by various workers 
p^ . However, in these works neutrino mass matrices are constructed using right-handed 
singlets along the line of the seesaw mechanism. 

In a previous paper |l^ we constructed a viable model of radiative neutrino masses 
with minimal SM matter in a five dimensional (5D) field theory on the orbifold S'i/(Z2xZ2) 
with bulk SU{3)w gauge symmetry. The SU{3)w unifies the SU{2) x U{1) electroweak 
symmetry of the SM |14|. The crucial observation is that the SM lepton doublet and 
charged lepton singlet can naturally be embedded in the fundamental representation of 
SU{?>)w- With Higgs fields in the 6 and 3 representations lepton number violating inter- 
actions can be constructed. Tree level masses are forbidden by orbifold projections and 
thus avoiding the fining tuning of VEV for small neutrino masses. Neutrino masses can be 
generated at 1-loop level. The scale of neutrino masses is small compared to the weak scale 
due to three factors: (i) the loop factor, (ii) the inverse of the compactification radius, R, 
which controls the volume dilution factor and (iii) small Yukawa couplings of the charged 
leptons. We found neutrino masses of order 0.01 eV without much fine tuning. Interest- 
ingly the solutions that satisfy the observed mixing parameters are found for the inverted 
mass hierarchy type if Yukawa couplings are not fine tuned. 

We continue this investigation for the GUT theory of SU{5) which is theoretically well 
motivated. This poses new challenges since the fundamental representation 5 unifies the 
lepton doublet e)^ with the right-handed down type quark dji. We found two options 
for the Higgs fields that can give rise to lepton number violating interactions essential for 
our mechanism. Strictly speaking the concept of lepton number is not fundamental in these 
theories; however, we find it convenient to use it both for guiding model construction as 
well as in navigating the tight constraints imposed by the many rare decay experiments. 
The Higgs fields are the 15 or 10 representation. Another problem we encountered in this 
construction is to maintain gauge coupling unification since it is well known that exotic 
Higgs contributes to the running of the coupling constants. This solved by introducing 
bulk fermion fields or additional Higgs fields. 

This paper is organized as follows: In section 2, we will review the setup of 5D orbifold 
SU{5) theory. The details of the model are given here. The construction of the Majorana 
neutrino masses through the Yukawa interaction of bulk 10 or 15 Higgs is made explicit 
in section 3. The gauge unification question will be addressed in section 4. The rich exotic 
Higgs sector leads to interesting phenomenology and will be discussed in section 5. Finally 
we give our conclusions in section 6. 



2. 5D Orbifold SU{5) Model 

We begin by a brief review of the 5D SU{5) unification model defined on the orbifold 
51/(^2X^2) |l^ with coordinates x^(;U = 0, 1, 2, 3) and the extra spatial dimension which 
is denoted by y. The circle 5i of radius R, or y = [— vri?, vri?], is orbifolded by parity 
P : y ^ —y transformation. The resulting space is divided by a second acting on 
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y' = y — ttR/2 as P' : y' ^ —y' to give the final geometry. On this orbifold, the Fourier 
decomposition is summarized in Table 1. 

Under the Z2 x Z2 transformation, there are two fixed points at y = and y = itR/2 
denoted by ys and ya respectively. The following two parities matrices are chosen for the 
parity transformations: 

P = diag{+ + + + +}, 

P' = diag{ + +}. (2.1) 

These determine the Kaluza-Klcin decompositions of a generic bulk field A(x, y) (sec Table 
1). The parity assignments of a bulk field arc chosen by phenomenological considerations. 
Indeed for a given multiplet different components can have different parities under Z2 x 
In particular, for the SU{5) gauge fields Am{M = /x, y), the following parities are used 

Z2:A^^ +PA^p-\ Ay ^ -PAyP-^ 

Z'^-.A^^ +P'Ai,P'-\ Ay ^ ±P'AyP'-\ (2.2) 

where we have suppressed the group index. Since the fifth components vanish at the ys 
fixed point they will not enter the low energy effective theory. With this parity assignment, 
when decomposed into the SM subgroup G = SU{3) x SU{2) x f/(l), the SU{5) gauge 
bosons have following components: 

24 = (8, 1,0)++ + (1,3,0)++ + (1,1,0)++ + (^3,2,-0 ^^'^^ 

where hypcrchargc is normalized to Y = Q — T3. The parities {P,P') are shown as sub- 
scripts. The first three terms are the SM gauge bosons. With the assigned (++) parities 
they will have zero modes. On the other hand, the last two terms represent the gauge 
boson Xs and Ys which have no zero modes since their parities are (H — ). They are KK 
excitations. At the symmetry broken brane, ya, even these Kaluza-Klein(KK) excitations 
vanish. So at ya, the symmetry is always G. On the other hand, at the symmetric brane, 
ys, SU{5) is still a good symmetry. At the low energy, only the zero modes of G can be 
observed. When energy is higher then 1/R, the appearance of KK excitations of SU{5)/G 
plus the KK excitations of G will restore the full SU{5) symmetry. By geometrical con- 
struction, the gauge X, Y gauge bosons are massive since they are KK excitations and the 
lightest ones have masses ~ 1/i?. At the bulk SU (5) symmetry is broken by the orbifold 



(P. P') form iiiMSs 

^++) ^[^o(x) + V2E„=i^^„+(^)cos^] 2n/R 

;^[^Sn=i ^r-i(^) cos ^] (2n - 1)1 R 

(-+) ^W^T.n=i ^2-nti(x) Sin ^^^] (2n - 1)/R 

Table 1: KK decomposition of a bulk field A{x, y) with parities (P, P'). 
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Figure 1: The proton decays through exchanging the X and Y gauge bosons. The subscripts 
denote SU{5) representations. 



boundary conditions. So there is no need to introduce the 24 as in the 4D case. Thus, the 
triplet-doublet problem |jl^ is solved naturally To break G to SU{3) x U{1) and to 
generate fermions masses, one still needs one Higgs in SU{5) fundamental representation. 

The SM fermions can be grouped into the standard 5 and 10 ( henceforth denoted as 
F and T respectively ) representations: Fi = (DfjLi) and Tj = {U^,Qi, Ef) where we use 
left-handed chiral fields with obvious notations and i = 1,2,3 is the family index whereas 
the superscript c denotes charge conjugation. In a 5D theory they can either be bulk fields 
or brane fields. If the matter fields are placed at the symmetric brane they will enjoy 
the merits and suffer the drawbacks of SU{5) GUT: the quantization of hypercharge, the 
unification of gauge couplings, the prediction of mass ratio of down-type quark and charged 



lepton, and the baryon number violating decays |17]. In a 5D model the possibility that 
some fermions can be bulk fields and plus the existence of different fixed points to locate 
fermions can be used to retain the successful mass relationships and avoid some difficulties 
such as the proton decay problem. For simplicity we shall assume that the Fi are localized 
on a brane at ys- Hence, there will be no KK excitations of neutrinos. 

Because the prediction of mfj/nir agrees quite well with experiment it is reasonable to 
localize both T3 and F^ on the brane at ys- On the other hand, we can make Ti a bulk field 
so as to avoid rapid proton decay, see Fig.|. This is so because both Ti and the X, Y bosons 
are bulk fields and their interaction must honor the conservation of KK number; i.e. the 
X, Y bosons are KK modes they do not couple to two up-quark zero modes. Therefore, all 
the sub-diagrams vanish in the 5D model. Similar diagrams which lead to baryon number 
violating neutron decays are also forbidden by the KK number conservation and hence 
matter stability is not violated. 



As pointed out in |18] the parities of bulk fermion fields pose further complications. If 
Ti has parities (++), then the parities of its components will be (f/f (++), Qi{-\ — ), Ef{++)). 
Clearly the quark doublets will not have zero modes and hence not acceptable. To complete 
the SM particle spectrum another T[ of parity (+— ) to be introduced such that the combi- 



-4- 



nation of the zero modes of (C/f , Q'l, E^) and brane field Fi constitutes tlie first SM family. 
In doing so the SU{5) mass relation will not be obeyed which is good phenomenologically. 
Henceforth we shall adopt the proposal that both Ti and T2 are bulk fields and no SU (5) 
mass relations holds for the first two generations. 

Letting Ti and T2 go into bulk immediately implies a hierarchy between 4D effective 
mass of the third and the first two generations. The 5D Yukawa interaction can be written 
as 

with group indices suppressed and M* denotes the fundamental scale. Because Ti and 
T2 live in the bulk, by naive dimension analysis, the 4D effective Yukawa coupling for 
down-type quarks naturally exhibit the following hierarchy patterns 



Ud ~ Cd 



(eel 

eel \ (2.5) 



where the volume dilution factor e ~ {Mc/M^,)^/'^ ~ 0.1 and q is a common factor. For the 
up-quarks, we adopt the following structure for 4D effective Yukawa couplings as proposed 
by 0: 

/e2 £2 e 



Vu 



e^ e^ e \. (2.6) 
e e 1 



The volume dilution mechanism still works for the Yukawa interaction between one brane 
and one bulk fermions, i.e. yui3,yu23,yu3i and 2/^32- But one needs to assume the Yukawa 
couplings between Ti,T2 to be two order of magnitude smaller then the others in order 
to have the structure of Eq. ( |2.6D . Given the totally different nature of Ti , T2 from other 
fermions, this assumption is an ad hoc one and we have no deeper understanding why this 
has to be so. 

The 5D bulk lagrangian is Z2 x symmetric. The question now arises: what parities 
should be assigned to the brane fermions? In our case the brane fermions will reside at 
y = 0, or equivalently y = nR, brane. They have definite parities under Z2 transformation, 
such as 

Z2:^-5^ +P^5, ^-10 ^ +P^^ioP (2.7) 
: ^-5 ^ iP'^-g, ^-10 ^ ±P'^^ioP' (2.8) 

So there are two ways of how the Z2 parities can be assigned to brane fermions: (1) 
(T, F) = (Q, U, E, D, L) = ±(+--, -+) and (2) (T, F) = (Q, U, E, D, L) = ±(+ — , +-). 
In order that the 5D lagrangian shall respect the symmetry of the theory we can simply 
make the follow manipulation on the brane term: 



j d'x I dy^{6iy) ± 5{y - 7rR)}L±, (2. 
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where subscript + and — stand for Zg-even and -odd respectively. Since our effective 
physical space is restrict to y G [0, 7ri?/2], the choice of parity for the brane fermions do 
not make any difference. 

Now we turn to the discussion of electroweak symmetry breaking. This is done by 
the bulk scalar in the 5 representation, H^, to reduce G to SU{3) x U{1)q and to give 
charged fermion masses as seen in Eq.(p.4|). However, contains the color Higgs He 
with quantum number (3,1,-1/3) and the ordinary SM Higgs doublet Hw The (H — ) 
parity of forbids the existence of zero mode and hence no light color scalars on the 
ys brane. Denoting = {h^ ,hP / \/2)'^ as in SM, the Higgs doublet develops VEV 
{H'^) = (0, v^^'^ / \/2) which breaks the electroweak symmetry as usual. Integrating out the 
fifth dimension, the effective gauge coupling can be used for relating the 5D and 4D gauge 
couplings, g4 = g/\/ ttRM* . The resulting W boson mass is found to be = g2TTRvf/2, 
or (27ri?z;^)^/^ = tjq ~ 250 GeV. With the above parameter substitution, we arrived at the 
effective 4D Yukawa interaction relevant to charged fermions masses: 

= yu,ij-^u'RiULj + yd,ij^{d' Rid'^j + e'Rje'ii) + H.c. (2.10) 

= u'rMuu'l + d'nMdd'L + e' RMee'^ + H.c. (2.11) 
Where the mass matrices can be identified as: 

Mu = Md = M'^,= ^Vd. (2.12) 

The mass matrices are diagonalized by bi-unitary rotation: 

UrMuUI = M'':'^, VRMdVl=Mf, (2.13) 
VlMeVS=Mf3, Vckm = uIVl (2.14) 

and the leptons can be rotated into their mass eigenstates by 

L, = iV^^jL'^, em = {Vl)^,e'R^. (2.15) 

Note that the neutrino is massless at this point, so we have the freedom to apply left-handed 
rotation to the entire lepton SU{2) doublet and make the charged current interaction 
diagonal. Also the Eq.(2.5,2.6) can not be taken literally otherwise the first two generations 



are degenerated. It shall be understood that every entity exhibits a small correction to the 
leading term exhibited. Numerical check shows that by allowing a 10% correction of each 
entry, it's easy to get the resulting mass hierarchy: 

m,:ms:md = mr:m^:me^l:e:e\ (2.16) 
mt : TTT-c : m-n ~ 1 : : e*^, {Vus, Vet, Kb) ~ (e, e, e^)- (2.17) 

In our numerical experiments about one third of the statistical samples provide CKM 
matrix very close to experimental values. Hence, these two general mass matrix patterns 
successfully yield the observed charged fermions mass hierarchy and CKM mixing angles. 
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Since the lepton mixing is crucial in the mechanism for generating neutrino mass which 
to be discussed in section 3, more insight of the mixing matrix is required. It's easy to see 
that Vl ~ diag{l, 1, 1), so the CKM is mainly controhed by the mixing of up-quark, Ul- 
The Vr is to diagonalize the almost uniform mass square matrix 




(2.18) 



For the leading approximation, the Vr is 

Vr 



V2 1 1 

V3 V(i Vti 
111 

\ V3 V3 V3 / 



(2.19) 



and it is equivalent to all its variants with any column permutation or arbitrary row sign 
flips because 

MdMl = SiMdMlSi, {Mff = G^{MffGi (2.20) 



where 



Si 



(o 


l\ 




/o 1 o\ 




^1 o\ 





1 


,52 = 
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,53 = 
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00 J 




[ool) 




[oioj 



(2.21) 



correspond to the interchanging of Vrs 1-3, 1-2, and 2-3 columns. And Gi = diag{—l, 1, 1), 
Q2 = diag{l, —1, 1) and ^3 = diag{l, 1, —1) are to flip sign of first, second and third row 
of Vr. This discrete symmetry is broken when one allows for small corrections to A^^ as 
discussed before. Sometimes the first two indices switch to get the right order of charged 
lepton mass, rUe < m^. 

From our discussions above, we see there is a 5D non-SUSY orbifold SU{5) model setup 
which can give a good charged fermion mass hierarchy and mixing angles. Furthermore, 
this framework preserves the advantages of SU{5) GUT and avoid the major obstacles in 
the traditional 4D theories. Since it's not SUSY, the unification of coupling constants is 
nontrivial. We will defer the discussion of RG running and gauge unification to section 4. 
Now we turn to the neutrino masses problem. 



3. Neutrino Masses 

It is well known that the symmetric 15 Higgs can be used to generate neutrino Majorana 
masses and break {B — L) in conventional SU{5) theories |1£]. For the orbifold version, we 
closely follow the study of [|13| where a model of neutrino masses in 5D orbifold SU{3)w 
unified theory without right handed singlet was constructed and the symmetrical 6 scalar 
filed was responsible for the neutrino masses. In the present case because 5 x 5 = 10 + 15, 
besides the 15 Higgs, the anti-symmetrical 10 Higgs can also work for generating the 
Majorana masses through one loop diagram, see Fig.0. 
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Figure 2: The 1-loop neutrino mass through {&)H^ Hi^Ht^ and {h)H^^ HiqH^ couphng. 



Two things differ from the SU{3)uj case studied by Firstly, in the SU{3)uj theory 
both 3 and 6 contribute one Higgs doublet. One of the linear combination is the would be 
Goldstone boson to be eaten by the SM gauge bosons and the orthogonal linear combination 
is the required physical charged scalar that appears in the loop. In SU (5) case , the 10 
and 15 break into 

15 = (6, 1, -2/3) + (3, 2, 1/6) + (1, 3, 1) (3.1) 

Pl5 Ci5 T 

and 

10 = (3, 1, -2/3)^ + (3, 2^1/6)^+ (1, 1, ly (3.2) 

Pio Cu) S 

in terms of SM quantum numbers. For notational simplicity, the symbols P, C, S, T have 
been introduced to indicate the specific component of 10 or 15 as shown in Eqs.( |3.l| , ^.l| ). 
We see that neither 10 nor 15 contain a Higgs doublet component. So in the SU{5) case 
we have to introduce another 5' Higgs for the radiative mechanism to work. 

Secondly, due to the SU{5) symmetry which imposes strong constraints on model 
building the parity of exotic Higgs sector needs careful examination. Now we discuss the 
various possible parity assignments of 15 and 10. 

We first assign the 15 parity to be (H — ). The 15(H — ) will decompose into -Pi5(H — ) + 
Ci5(++) + — ). By doing so, the triplet component T has no zero mode and hence 
no VEV can be developed. Hence, there is no tree level neutrino Majorana mass. Also, 
the hierarchy of the VEVs of SM Higgs doublet and the T Higgs triplet required by the p 
parameter measurement will be naturally avoided. In this case, the C15 color Higgs will 
improve the unification of non-SUSY SU{5) [pO| . 

There are several shortcomings of this parity assignment. The Higgs boson running in 
the loop (see Fig.^) are both KK modes due to KK number conservation. So the resulting 
neutrino masses are too small due to the suppression of KK mass 1 /R. Furthermore the 
parity of the extra 5' Higgs has to be (H — ) such that the bulk triple Higgs interaction 

-^h'^H,,„,HI+H.c. (3.3) 

is Z2 invariant. Due to the above Higgs potential proton decay now will be induced via 
the mixing of the zero modes of C15 and the color component (3, 1, —1/3) in 5', see Fig.^. 

To satisfy both the constraints of proton decay and the requirement that the resulting 
1-loop neutrino mass is of oder .01 eV, the zero mode mass of C15 need to be as high as 
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Figure 3: The proton and neutron decays through the mixmg of (3, 1, —1/3) color Higgs in 5' and 
the (3, 2, 1/6) color Higgs boson in 10 or 15. 

10^^ GeV and the triple Higgs coupling strength m is about 10^^ GeV. This scenario is 
disfavored by the strong interacting scalar sector and a very heavy scalar. 

Next, the parity of 15 can be chosen as ( — 1-). Then the decomposition is Pi5(— +) + 

Ci5( ) + r( — h). Again, the triplet T has no zero mode, so naturally avoid the VEV 

problem and the tree level neutrino mass. Now parity of the extra 5' Higgs has to be (— +) 
correspondingly. But the lepton number breaking component vanish at y = brane, so 
it doesn't couple to the 5 living there. For this scenario to work, the 5 fermions have to 
be placed at y = 7ri?/2 and enjoy no SU{5) symmetry. Unfortunately, even doing so, the 
resulting neutrino mass pattern is not right. 

The third parity assignment is ( ). This leaves the T components totally vanish at 

both branes. It does not work in our scenario where the Fi are localized on the brane. 

So we are only left with the (++) assignment. Explicitly, it satisfies the following 
transformation properties 

PHi^P-^ = +Hi5, P'Hi^P'-^ = +Hi^. (3.4) 

Now there are ^15^^^++), {C^i''\C^i'\+-) and (r+2, T+i, rO)^(++) Higgs. The 
couples two neutrinos which could give them the tree-level Majorana mass. To avoid 
that, we assume that it is a regular scalar and does not develop a VEV. Also, the zero mode 
of Pi5 must be heavy to suppress the K — K mixing, which can be gleamed from Fig.^ 
This gives a satisfactory neutrino mass matrix and other phenomenological constraints. 

Most of above discussions on the parity assignments can be applied to 10 so we will 
not repeat the analysis here but note that 10 doesn't acquire VEV. The suitable 10 is of 
(++) parity under Z2 x Z'2 too. It has components Piq^^{++), {Ciq^^ ,C^q^^){-\ — ) and 
S+\++). 

Let's summarize the setup of the scalar sector: (1) One extra 5' Higgs. (2) Either 15 
or 10 scalar with parity (++) under {Z2 x Z2). (3) The scalar 15 must not have VEV 
to avoid excessive fine tuning. The zero modes of scalar 15 or 10 must be heavy to avoid 
large contribution to K — K mixing. 

The necessary lepton number violation terms arise from the brane Yukawa interactions 



C-Yib = 5 (y) 



fl5 



V2M* 5j 15 



(3.5) 
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for 15 where h\^^^ = H^^"^^ and 



>Cyio = 5 (y) 



flO 



(3.6) 



for 10 and H^q^^ = —H^q^^. After integrating out the extra dimension, we then param- 
eterize the relevant 4D effective interaction as: 



C D -J^e'-iu'S+' + ^ e'-iv'T+^ + H.c. (3.7) 



in the leptons' flavor basis. It is easy to see that flj' = fjf and flj' = —fji- The entries 
of the neutrino Majorana mass matrix can be approximately calculated by ignoring the 
lepton masses in the propagators and are given as follows: 

where /io/i5 jg ^j^g Yukawa coupling of either 10 or 15, the mass of charged lepton-A: 
and Ml, M2 the masses of Higgs running in the loop. And /'^ is the effective 4D Yukawa 
coupling for the relevant physical charged Higgs interaction: 

f-^e'ji^,^',H- +H.C. (3.9) 

It is more convenient to express the neutrino mass matrix in the basis of charged 
lepton's mass eigenstates. Alone the line of the discussions in Sec. 2, we expect the lepton 
Yukawa coupling (/ ^)^ to exhibit the same pattern as yd in Eq. (|2.5D . As we rotating 
the charged leptons into their mass eigenstate, the Yukawa matrix / ^ is more or less 
diagonalized in the same time and it is proportional to the mass of charged lepton: 

'5 X 92mi 

In the lepton mass eigenbasis, the Yukawa couplings for 10 or 15 can be obtained by the 
following transformations 



Only Vfj is involved since we have invoked the SU{5) mass relation between the d'^ and e. 
Note that the antisymmetry of flavor indices i,j of flj^ and the symmetricalness of flj" is 
not affected by these transformations. 

Therefore the Eq.(|3.8|) is further simplified to 



2 ^10/15 , ^9 

^^^^^ ^ 6A.^M^inRM*)Mf-Mi M? + ^ ^'''^^ 



Thus, the matrix to diagonalize Eg. ( 3. 12 ) is the MNS matrix and we can directly use the 



phenomenological results of neutrino oscillation studies. With vq ~ 250 GeV, rrii ~ ruj 
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/10/15 ~ 1, Ml ~ 10^ GeV, M2 ~ 200 GeV and (ttRM*) ~ 100, the overall neutrino mass 
can be estimated to be 0.13 x (/^'^/^^)(m/TeV) eV. It can easily give neutrino overall mass 
~ 0.01 eV by adjusting m ~ 0.1 TeV such that the triple scalars is weakly interacting, see 
Eq.(O). 

For the case of 15, the neutrino mass matrix takes the form: 



Ml' ~ ^ 



1 

M 



( 



mi f 



f 



\ffl{mi + mi)/2 fM{ml + mi)/2 ,33, 

where M is a mass to be determined by Eq.( |3.12j ). Basically, it can yield any of the 
acceptable neutrino mass matrices. Here, we give few examples to illustrate the richness 
of resulting neutrino mass pattern by using 15. 
For instance, one solution is 

/ 0.4959 0.2012 -0.7524 \ 
= 0.2012 0.2717 -0.4316 (3.14) 
y -0.7524 -0.4316 1.2360 j 

which leads to a normal hierarchy mass pattern: 

/ 3 X 10-5 0.4233 -0.0271 \ 
ruo^M^ = 0.4233 1.0 0.9319 (3.15) 
Y -0.0271 0.9319 1.1493 / 

with 01 = 37.4°, $2 = 44.04°, S3 = -0.13, AQ/Aatm = 0.017. mo sets the over- 
all neutrino mass scale which is determined by factors in Eq.( 3.12| ). In this example, 
mo = 0.05(m/0.15TeV) eV. Notice that the first entry will predict an unobservable rate 
for neutrinoless double beta decay. We have searched numerically for larger values of 
and found solutions with first entry as big as ~ 0.1. But they require fine tuning 
of Yukawa of order 10~^. Even this value is below the sensitivity of the next generation of 
these experiments. 

For neutrino mass matrices of the inverted hierarchy type, this mechanism for gener- 
ating neutrino masses requires fine tuning of Yukawa couplings f'^j^ which we will discuss 
below. Following the classification in [0, ^, we explore the following two leading patterns 
which lead to inverted hierarchy: 



[m: 



ml)/2 



+ ml)/2 fll{ml + mi)/2 
fll{ml + m^)/2 



2 

e 

15 
22 

2 I ^2 



,m. 



fl^ml 



(3.13) 



mo 
V2 




M^"^ ~ mo X 



/l 0' 

OH 
VoH, 



(3.16) 



Both of the leading patterns give bi-maximal mixing angles which confiict with the current 
experiments. Some small entries in the structure zeros as well as small perturbations to 
the leading terms are necessary to accommodate the experimental data. 
Let's examine two examples: For B\ type, one solution for My \s 

/ 0.42 1 0.922 \ 
Mu = 1 0.097 -0.464 (3.17) 
»^ 0.922 -0.464 0.006 j 



V2 
mo 
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which gives 6*1 = 36.57°, 6*2 = 42.43°, S3 = 0.06 and Aq/A 



atm 



0.021 and imphes 



JBl 



/ 3 X lO"*^ 7 X 10-5 
7 X 10-5 0.6667 
y-5 X 10-5 -0.4715 



-5 X 10-5 
-0.4715 
0.3335 



(3.18) 



One example of B2 type mass matrix is: 



mo 



( 1 0.02 
0.02 0.49 
-0.01 0.5 




which gives Qx = 34.36°, 62 = 45.29°, S3 
couphng in the lepton weak basis is 



0.021 and Aq/A 



atm 



(3.19) 



0.030. The Yukawa 



fl5 
JB2 



/ 6 X 10-6 
-3 X 10-6 
y-6 X 10-6 



-3 X 10-6 
0.6667 
-0.4714 



-6 X 10-' 
-0.4714 
0.3333 



(3.20) 



From the above examples we see that for both cases elements of the first row and first 
column are much smaller than the rest. This amounts to fine tuning of the model. We 
note in passing both inverted hierarchy types can be obtained from a Yukawa matrix with 
the leading structure that looks like 



15 




/ 



where x denotes a small number of order 10 5. 

For 10, the resulting neutrino mass matrix structure is: 

/ 



(3.21) 



1 

M 







V/: 



m 



m: 



mi)/2 
ml)/2 f^l{ml-ml)/2 



2 rtTO 2 

e 
2 





ml)/2 
ml)/2 



(3.22) 



Again, M is some mass to be determined by Eq. (|3.12 ). The diagonal zeros are the result 
of antisymmetry of the 10 representation. It can only give Bl type mass pattern, namely 
of the inverted hierarchy type, if the Yukawa exhibit the following hierarchy: /J{^2 ■ /lis • 
~ 1 : m'i/m'^ : m^/m^ which implies that 



?io 








1 


1 








1 





1 






V2 


v 


1 


1 







V2 



(3.23) 



The hierarchy is also a refiection of the approximate Lr + L^ — L^ symmetry. It's impossible 
to get realistic neutrino masses by 10 alone at this leading order. However, small diagonal 
entries can be generated through 2-loop correction which we shall leave for future works. 
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Another obvious mean is to include 15 as the perturbation source but this will lead to 
non-minimal models. 

Hence, it is interesting that 15 gives a normal hierarchy whereas the 10 give inverted 
hierarchy of neutrino masses without additional symmetry or fine tuning of their respective 
Yukawa couplings. 

4. Gauge Unification 

It is well known that the minimal 4D SU{5) unification is ruled out by current data. Even 



the supersymmetric version pl| , p^ j is disfavored |23]. With the extended Higgs sector 
required to generate neutrino masses we expect the unification of the gauge couplings will 
be a challenge. Since our model is a 5D one albeit non-supersymmetric, the renormalization 
group (RG) running of the SM gauge couplings crosses several thresholds and scales. We 



follow here the methods developed by the works of p5| , 26]. Below the compactification 
scale, 1/R, we have 4D eff^ective field theory with thresholds crossings from the 15 or 10 
mandated by our model of neutrino masses. This will constitute an intermediate scale 
between that of electroweak breaking and 1/R. Between 1/R and the unification scale one 
has power law running of the 5D gauge theory. Clearly the problem of gauge unification 
now becomes highly nontrivial and requires careful treatment which we detail below. We 
will only discuss RG running and unification up to one loop. 

The first task is to establish the lower limit on the mass of the 10 or the 15, generically 
called Mp. We found the most stringent limit for Mp comes from the Kl — Ks mass 
difference mediated by exchanging color Higgs components of -Picis, see Fig.^. 

d d d ^ s 

^15 Pi5,Pio 
(a) s s (6) s d 



Figure 4: The tree level contribution to Kaon mass difference by exchanging the (6, 1, —2/3) Higgs 
in 15 representation ((a),(b) ) or the (3, 1,-2/3) Higgs in 10 representation ( (b) ). 



The A 5 7^ relevant terms are: 



flO 
Jij 



al3] 



+ 



fl5 
Jij 



-J?dfp}f^ + H.c. 



(4.1) 



VirRM* 

where a, (5 are the color indices and it is understood the above expression shall be made 



SU{2>) invariant (for details see |27|). The Yukawa interaction can be rotated into down 
quarks' mass eigenstates again by Eq.( p.ll ) because the SU{5) symmetry relates left 
handed leptons and right handed down quarks. After applying the Fierz transformation, 
we arrive at the eff'ective AS = 2 operator from exchanging the P component of 15: 



AS=2 



1 1 



2M|, 



fl5 fl5t I I fl5|2 
Jll J 22 "T I Jl2 I 

ttRM* 



{dl^,sl){dirsi) + H.c. 



(4.2) 
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By using vacuum insertion approximation, the resulting kaon mass difference can be esti- 
mated to be: 

AMk ~ - fKmKRe ( lEl^±}Il!l] . (4.3) 

Tlie formula applies to 10 except that fl^ has no diagonal element. Plugging in the 
experimental values: fx ~ 0.16 GeV, mx = 0.4976 GeV and AM^ = 3.48 x 10"^^ GeV 
we obtain the lower limit for Mp: 

Mp > 1.11 XIOM ^11^22 +1/12 I I \ Qgv. (4.4) 

\^ nRM* j V3-48 X lO-iSQeVy ^ ' 

where AM^ is the mass difference due to color scalar contribution. Similar analysis can 
be naively extended to AMb by substituting Jk ^ Jb and Mk ^ Mb which gives a less 
constrained limit: 

^/;f/3r+i/;ih"V ami 



Evidently, the limit goes down as the absolute value of specific Yukawa couplings of 15 or 
10 is lower. 

If we take the above number as the intermediate scale the question arises whether 
these color Higgs will induce rapid proton decay. It can be checked that there are neither 
tree level nor one loop contributions. Hence, with a relatively low value of Mp ~ 10^ GeV, 
obtained by assuming that ttRM* ~ 100, AMx is saturated by the contribution of exotic 
scalars and the extreme case 

1/10/15 1 ~ 1, will not run afoul of proton stability. 
The one loop gauge coupling RG running after passing various thresholds can be 
written as 

= ar^(Mz)- 4^ In -^-^ In ^ 



2tt Mz 2tt Mh 

N° N'^ 

_ of fiR 
27r ^ ■ 

71=1 n=l 



-^Vln^^-^Vln^ (4.6) 
27r ^ 2n - 1 27r ^ 2n ^ ' 



The last two a terms are the KK modes contributions when energy scale crosses 1 /R. The 

integers N° and A^^ are the highest odd( (+— ) and (— +) ) and even( (++) and ( ) ) 

KK excitation level below the scale ^ 

2Are - 1 , 

"fT"^' R ^ ^ 

For the SM, the beta functions are well known: ap''' = (4, -10/3, -7)+ni^x(l/10, 1/6,0) 
where uh is the number of Higgs doublet zero modes. At scale Mz, deriving from 
a{Mz)-^ = 127.934(7) and Ow = 0.231113(15)[p8|], we have 



1 3 cos 9w , X 1 sin Qw , 

a7^ = = 59.031(35), = = 29.568(17) 4. 

5 a a 
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and as{~l) = 8.53. In our model, the zero modes of 10 or 15 Higgs has a large mass Mp 
imposed by phenomenology. Below Mp, only SM particles go into the beta function. As 
one passes the threshold, Mp < ij, < ^, the zero modes of 10 or 15 will contribute to the 
beta function. For 15, it has zero modes Pi5(6, 1, —2/3) and T(l, 3, 1) so 

(ai^4^ar)=f^,yV (4.9) 



^15' 3' 6, 

For 10, the zero modes are Pio(3, 1, —2/3) and S{1, 1, 1) hence 

«,4°,4°) = (^,o,i). (4.10) 

Finally, when the scale is over the compactification scale, various KK excitations gradually 
come in and contribute to RG running. For very high KK levels, summing over the KK 
tower roughly gives the power law running of coupling constants. This power law behavior 
can be easily understood by adding up the KK excitations level by level. By simple algebra 
and Stirling's approximation, the KK sums are 

n=l 

~ N° ln(/xi?) - AT" In 2Ar<' + - In + C>(-^) (4.11) 

n=l 



- N^'lnd^R) - Nnn2N'' + N"" -\nV2^ + 0{^). (4.12) 

When /i S> 1/-R, N° ~ A^'^ ~ fJ'R/'^, the contribution of all the KK excitations to a~^{fi) 
can be approximately calculated as follows: 

(4.13, 

where Oj = + a° is the sum of the beta function from even and odd KK components. As 
H being higher then the compactification scale the full SU{5) symmetry start to emerge. So 
it is not surprising that any complete SU{5) multiplet gives equal amount to three gauge 
running, ai = 0,2 = a^, which will not affect the gauge coupling unification. As unification 
is concerned, only the even (or equivalently the odd ) KK components matter. Explicitly, 
we list the beta functions of all even KK contents in Table 2. 

Note that the fifth components of the KK X, Y gauge fields come in as real scalars 
which give one half of the contribution of complex scalars. Also note that the combined 
result of bulk filed T and T' are same to three coupling running so it will not affect 
unification. The contributions of KK modes can be summarized as: 

,~e ~e ^1 32\ / 1 1 ^ 

(a„a„a3)=^-,--,--j+n5X^-,-,0 
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a 


Even Components 


a? 


z 


a 


-"5 


1 /fi 

i/D 


-"10(^1, 1/ 


1 /I n 


1 /6 
i/D 


n 

u 


TT 

-"15 




M5 (0' 1' -2/3) 


o/ io 


n 
U 


^ /a 
0/D 






i (1,3,1) 


o/o 


2/6 


U 




1/2 


A~o'^'(3, 1,-2/3) 


4/15 





1/6 






5(1,1,1) 


1/5 












0^^(8,1,0) 
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(1,3,0) 





-22/3 





bU (5) Gauge 


—35/2 


^^^(1,1,0) 















X5,y5(3,2,-5/6) 


5/12 


1/4 


1/6 






X5,y5(3^2,5/6) 


5/12 


1/4 


1/6 


Tbulk 


2 


[7^(3,1,-2/3) 


16/15 





2/3 






^^(1,1,1) 


4/5 








rpl 

^bulk 


2 


Q'l(3,2,1/6) 


2/15 


2 


4/3 



Table 2: RG contribution from various KK excitations. 



and 

35 1 7 1 , , 

a = — - + -ns + -ni5 + -nio + 4nT (4.15) 
zoo 2 

where ns, nio, nis are the numbers of bulk 5, 10, 15 Higgs and tit is the number of 
generation of ten-plet fermion in bulks respectively. With the minimum particle content 
we cannot obtain gauge unification as can be seen explicitly in Fig. 5. 

To make this model work for unification, we propose to use vector fermions. The 
vector nature is to ensure no new anomaly is induced. They can be either ng pairs of 
heavy quark doublets, Qh + Qhj coming from of Th{-\ — ) + Th{-\ — ) or pairs of heavy 
lepton doublets Lh + Lh from bulk til pairs of Fh{++) + Fh{++). The new Fh + Fh 
must be engineered not to mix with the -^1,2,3, otherwise the neutrino mass matrix will be 
spoiled. For this reason we favored heavy quark pairs. 

For simplicity, we assume they all have a common zero mode mass, denoted as Mql, 
which can be adjusted to achieve unification. When Mql < fx < 1/R, the threshold effect 
of Qhs 

(Aai, Aa2, Aag) = 2nQ x 1, ^ (4.16) 
should be considered and the extra term 



should be included in Eq.(4^). As the scale larger then 1/R, their KK modes give extra 



contribution to d'^ and d° as 

Aaf = 2Aai, (Aa?, Aa^, A5°) = Ang x (^^, 0, ^ ) (4.18) 
where another factor 2 is due to that 5D fermions are vector like. 
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It turns out that a simple solution is with ng = 1. Now we have two 5 Higgs, one 
15 Higgs with mass Mp = 10^ GeV to evade the constraint from AMk and one pair of 
Qh + Qh- If ^6 require that (ttRMgut) = 100, then the solution can be found to be 
Mql = 18.03 TeV and Mqut = 8.94 x lO^^GeV, see Fig.6. 

Note the splitting of even and odd KK contribution slows down the convergence and 
push the unification scale to higher end. It is interesting to note that the volume factor, 
~ (ttRMgut), is related to Mql. A higher Mql provides smaller volume factor and 
lower Mgut- For example, to have (ttRMgut) = 1000 it requires Mql = 6.0 TeV and 
Mgut = 2.97 X 10^'^ GeV. We have the approximate upper bound Mql < 117.1 TeV and 
Mgut > 1-16 x 10^^ GeV from the consistency requirement that n > 1/R. 



Mp l/R Mq Mp l/R 




5 10 15 , 5 10 15 , 

Figure 5: The RG running with Mp = 10^ Figure 6: The RG running with Mp = 10^ 
GeV and l/R = 2.81 x lO^'' GeV. GeV, Mq = 18.03 TeV and 1/i? = 2.81 x lO" 

GeV. They converge at 8.94 x lO^^Qgy. 



We note in passing that unification can still be made by adding many bulk 5 Higgs 
with parity (++) if heavy fermion doublets are not used. However, this solution requires 
a large number(> 6) of extra scalars needed. We deem this to be an unpalatable solution. 



5. Phenomenology 

The neutron anti-neutron oscillation may be induced by the mixing of Higgs, see Fig. |^. 
The relevant 6-quark operator can be expressed as 

G^B=2u''u''d''d''d''d'' + H.C. (5.1) 

the parameter G/\b=2 is of mass dimension minus five. It can be estimated to be: 

^^^=^ {l/RYMli^RM*) < (^i?M^)^'^ (^-^^ 

in arriving the result, the mass Mp ~ 10^ GeV from the constraint of AMk, /s ~ /is ~ 1, 
and m ~ 1/i? ~ 10^^ GeV have been plugged in. It is safely within the experimental 
limit 

r^v-iv > 0-86 X lO^sec (5.3) 
or equivalent, G^b=2 < 3 x lO^^^GeV"^. 
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Figure 7: The Feynman diagrams for N—N oscillation. 

e e 



T 



±2 



±2 



(b) 



Figure 8: The Feynman diagrams for (a)muonium-antimuonium and (b) fi 3e transition by 
exchange Higgs T^^. 



The lepton flavor violating muonium {fJ.^e~ = M)-antimuonium (/i^e+ = M) tran- 
sition and rare decay /i ^ 3e can be induced at tree level by exchanging the T^^ scalars 
which belong to 15 Higgs, see Fig.|8|. The interaction of T with charged leptons is given by 



fl5 
Jij 



-.l^iljT^'^ + H.c. 



(5.4) 



^/■nRM* 

where Yukawa coupling /^^ is in the charged lepton's mass basis. Assuming there is no 
external electromagnetic fields and the mass difference of the mixed state, 6 , is small. The 
possibility of observing a the transformation of muonium into antimuonium can be written 
as P{M) ~ 5^/2r^ where T ^ is the muon decay rate. And the mixing can be estimated to 
be 

o ^15 

'22 



5 EE 2{M\'H^jj^\M) 
where a is the Bohr radius (a^^ 



9 f 15 f isf 



arue 



-Ka?M^{-KRM*) 
P{M) can be expressed as 

4 



P(M) . 1.1 X 10-^|§g (^) 



(5.5) 



(5.6) 



which is safely within the current experimental limit P{M) < 8.3 x 10 
The /Li — > 3e transitions can be described by an effective lagrangian 

^^11 -^12 j (e-c^^)(e-ec)+//.c. 



M^iirRM*) 



which leads to the branching ratio of ^ ^ 3e: 



Br{fi 3e) 



g^{7rRM 



*\2 



Mt J 



(5.7) 



(5.8) 
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The above equation can be expressed in term of AMx to eliminate the uncertainty of the 
absolute value of /^^: 

2 



Br{ij ^ 3e) = 2 



SM^AruK 



ArriK 



Mp 



I fl5t fl5 



12 I 



15t fl5 
J 22 



11 



|;15|2 



'12 I 



(5.9) 



where Am^ is the contribution to kaon mass difference by exchanging P scalar. We give 
explicit dependence of the ratio of Mp to Mt because we expect they will split after 
quantum corrections are taken into account. In the above we have used the relations 
between Yukawa couplings and the elements of the neutrino mass matrix, see Eq.( 3.13| ). 
So the ratio of Yukawa couplings can be replaced by the ratio of the corresponding elements 
in Mu- 



Br[^i 3e) ~ 3.02 x 10"^^ 



Am^ 
ArriK 



Mp 
Mt 



2miimi2 



(5.10) 



^7711177122 + (2^ 

It's straightforward to extend the analysis to r ^ 3/ transitions. Assuming that the 
hierarchy of the elements of neutrino mass matrix is smaller than factor 100, this model 
predicts 



Br{fj, — > 3e) : Br{T 3e) : Br{T 3fi) : Br{T — > /iee) : Br{T e^^) 



m 



12 



rrin 



m 



13 



77ln 



rrif. 



777: 



23 



/ 777, 



: - 



771t 



777, 



23 



777o 



rrif. 



m 



12 



(5.11) 



(-22 \iii-T/ rrii^ yrrir / 77732 \''^t^J ?7^ii?7722 
The extra suppression of mass ratio to the fourth power makes it very difficult to find 
experimental signal in r ^ 3/ decays. Hence, ^ 3e is the best probe of the model. 
But this model exhibits an interesting characteristic: only the neutrino mass matrices of 
Bl type have the chance to be benefited from large enhancement (77712/77722)^ such that 
/i — > 3e can be observed in near future experiments. 

The rare decays ^ 67, r — > /i7, b ^ sj etc can be induced by one loop diagrams 
which involve corresponding components of either 15 or 10. Clearly, these rare decays 
are useful tools for probing flavor physics. But, again, these process are suppressed by 
{M\Y / Mp)"^ plus the loop factor suppression which make their rates too small to be tested 
in foreseeable experiments. 

What about seeing the effects of the heavy quarks Qh + Since we expect their 
masses to be heavier then 10 TeV direct production of these particles is not likely at the 
LHC. One can look for their virual effects. Because they are vector like so they have no 



leading order contribution to the S and T parameters] 31]. 



6. Conclusions 

We have explicitly constructed viable models of neutrino mass matrix involving only three 
active SM Weyl neutrinos without introduction of singlet right-handed fermion in 5D orb- 
ifold SU{5) unification models. The crucial source of lepton number violation is the 15 
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or 10 bulk Higgs fields. This model preserves the orbifold solution to the triplet-doublet 
problem and avoids rapid proton decay mediated by leptoquarks at the tree level due to 
KK number conservation at the vertex coupling two 10 fermions which are bulk fields. 
However, at the one loop level proton decay can occur. Our estimate puts it safely within 
the experimental bound. 

In the class of models we studied the overall scale of neutrino mass is partly controlled 
by the Yukawa couplings of the exotic Higgs bosons to the brane leptons and triple scalar 
coupling and the mass scale of the exotic Higgs fields. It is of order 10~^ eV after all 
the phenomenological constraints are satisfied. We found that the 15 prefers a normal 
hierarchy; whereas the 10 gives an inverted hierarchy when two loop effects are taken 
into account. Otherwise only bi-maximal mixing can be obtained. Since the model is 
restricted to three Weyl neutrinos the mass matrix obtained is necessary Major ana. This 
is a necessary but not sufficient condition to induce neutrinoless double beta decay of nuclei. 
To do so the (11) entry of the mass matrix must be non vanishing. Generically without fine 
tuning of the Yukawa couplings, it is typically of order 10~^ eV which puts it outside the 
range of detectability in the next round of experiments. This is a reflection of the fact that 
normal hierarchy or Zee-like mass matrices are preferred. Stating this differently, if a value 
of {nieel ~ 0.1 eV is extracted from the next generation of experiment which we imagine 
to give positive signatures then we would require a 4 orders of magnitude hierarchy in the 
Yukawa couplings of 15 for this mechanism to work. 

As expected the introduction of exotic Higgs exacerbates the problem of gauge unifica- 
tion of SU{5). We found a solution to this in the 5D model by using a pair of vector bulk 
fermions 10 and 10. The masses of these fermions are in the 10 TeV range and unification 
occurs at ~ 10^^ GeV and the compactification radius 1/R ~ 10^^ GeV. On the other hand, 
these exotic scalars are promising source to give the universe enough matter anti-matter 
asymmetry. Detailed analysis of this will appear elsewhere [^]. 

In our study we have not considered in detail the origin of the charged fermion mass 
hierarchy. It is sufficient for us to put the third family on the SU{5) brane and the first two 
families 10 fermions as bulk fields. Perhaps this can solved by the split fermion scenario 
similar to |27] but this we leave for future considerations. 

It is well known challenge to test the physics of various models for neutrino masses. 
Currently the data do not distinguish the different classes of models let alone the many 
within one framework. For the seesaw model, generally the right-handed Majorana neutri- 
nos are of order of GUT scale and direct probe is out of the question. For the case of bulk 
right-handed neutrinos studies thus far done have indicated phenomenological tests are 
only possible if the compactification is very large 33 1. Since the bulk neutrino behaves 
like a sterile ones more structures in the oscillation pattern will be a good signature p^ . 
Currently no such structure is seen and more precise measurements will be needed. In our 
radiative scenario it is possible to see the effects of lepton number violating Higgs in rare 
decay. We found that once the mass Mp satisfies the constraint from AMk then most of 
the lepton and baryon flavor violating transition induced by 15 are far below the present 
experimental limits. This also puts it out of the range of direct detection in the near future. 
However, the decay /i — > 3e can be just below the current experimental limit. If found may 
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indicate the neutrino mass matrix is of Bl type inverted hierarchy which leads to a very 
low neutrinoless double beta decay rate. 
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